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Abstract 

Entanglement and quantum discord for a pair of nuclear spins s = 1/2 in a nanopore filled with a 
gas of spin-carrying molecules (atoms) are studied. The correlation functions describing dynamics 
of dipolar coupled spins in a nanopore are found. The dependence of spin-pair entanglement on 
the temperature and the number of spins is obtained from the reduced density matrix, which is 
centrosymmetric (CS). An analytic expression for the concurrence is obtained for an arbitrary CS 
density matrix. It is shown that the quantum discord as a measure of quantum correlations attains 
a significant value at low temperatures. It is shown also that the discord in the considered model 
has "flickering" character and disappears periodically in the course of the time evolution of the 
system. The geometric discord is studied for arbitrary 4x4 CS density matrices. 
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I. INTRODUCTION 



Quantum correlations in many-particle systems are responsible for the performance of the 
quantum devices (in particular, quantum computers) and give them significant advantages 
over their classical counterparts |2J. Until recently it was believed |3| that entanglement 
is responsible for quantum correlations |4j , and quantum devices can be created only on the 
basis of materials which can be easily prepared in entangled states [2]. However, it turned 



out that one can work out quantum algorithms 0, 0] which significantly outperform their 
classical counterparts by mixing separable (non-entangled) states. It also turned out that 
quantum non-locality can be observed in systems without entanglement [6j. From this, one 
can conclude that entanglement describes only a part of quantum correlations but not all of 
them. 

r-| 

According to the current understanding |l| , total correlations, quantum and classical, are 
determined by the mutual information [7] . The problem is how to divide the correlations into 
quantum and classical ones. The problem was solved independently by Henderson and Vedral 
8] as well as independently Ollivier and Zurek (9)]. The classical correlation in a binary 
system is determined by the total set of projective measurements carried out only for its one 
subsystem js|. Then the quantum correlation measure (the quantum discord) is determined 
as the difference between the mutual information and its classical part optimized over all 
possible projective measurements js, 9]. The quantum discord is determined completely by 
quantum properties of the system and equals zero for classical systems. 

A calculation of the quantum discord is a complex problem because it is necessary to 



execute many-parametric optimization. There: 
possible only in simple two-qubit systems 10, 



ore, analytic calculations for the discord are 



A special attention was devoted to a connection of the quantum discord and physical 



parameters of the system [3, |l3|. Such investigations open a direct way to experimental 
measurements of the discord. In the present article we investigate theoretically the entan- 
glement and quantum discord in a thin silicon film containing close nanopores filled with a a 
gas of spin-carrying molecules (atoms) with spin s = 1/2 in a strong external magnetic fielc . 
NMR line shape in such films was studied experimentally and theoretically 15, 16 ]. 
Fast molecular motion in nanopores do not average the dipole-dipole interactions (DDI) 
completely (space confinement) and the residual DDI is described by one coupling constant 



2 



which is the same for all pairs of interacting spins Il5| . Reducing the density matrix 
(which describes the time evolution of the system with the residual DDI) over all spins ex- 
cept the chosen pair one can obtain information about the pair entanglement and quantum 
discord in the system. The main goal of this article is the investigation of entanglement and 
quantum discord of spin pair in the nanopore filled with a gas of spin-carrying molecules 
(atoms). 

It is also significant that the reduced density matrix in the considered model is centrosym- 
metric (CS) [17]. In the present article we obtain an analytic expression for the concurrence 



181 ] of an arbitrary 4x4 CS density matrix. The geometric discord for an arbitrary 4x4 



CS-matrix is also calculated. 



II. THE REDUCED DENSITY MATRIX OF THE SPIN PAIR IN A NANOPORE. 

We consider a system consisting of iV spin-carrying molecules (atoms) of a gas with spin 
s = 1/2 in a closed nanopore in a strong external magnetic field B (Fig. [TJ). Molecular 
motions lead to a partial avera ging of the secular DDI, Hdz Q, [3, and the Hamiltonian of 



the residual (averaged) DDI is [15] 



ff* = y(3£-A (i) 



where D is the coupling constant, which is the same for all spin pairs [14J, [15], I 2 is the 
square of the total angular momentum, I a = Y^i=i If, If is the operator of the momentum 
projection of spin % (i = 1,2,..., N) on the axis a (a = x, y, z). In the initial moment of time 
the resonance 7r/2-pulse, turning spins by the angle 7r/2 about the axis y (Fig. [TJ, acts on 
the system. As a result, the equilibrium density matrix, po, in a strong external magnetic 
field is given by 

Po = |e^, (2) 

where Z = Tr{e /3/a: } = 2 7V cosh 7V | is the partition function, = is the inverse dimen- 
sionless temperature, ujq = jBq (7 is the gyromagnetic ratio) is the Larmour frequency, and 
T is the temperature of the system. The time evolution of the density matrix is described 



as 



p(t) = e- iHd ^poe iHd ^ = I e -^W a ^ 2 , (3) 
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Figure 1: A gas of spin-carrying molecules (atoms) in a nanopore in the external magnetic field. 



where we took into account that [I a ,I 2 ] = (a = x,y,z) and set a 



3D 



matrix ([3]) describes the free induction decay and NMR line shape 
our problems it is convenient to rewrite the density matrix p(t) as 19] 

3 



The density 



161 ] . In order to solve 



p(t) 



E 



a 



1,2,. ..,N 



where ^ = 0, 1, 2, 3 (j = 1, 2, N), J° = Ej (Ej is the 2 x 2 matrix 



i x i 2 

3 ' 3 



(4) 



7 y 

3 ' 



J| = //, and af i 2 ';^ {t) is a function of time. Since Tr(lf) = ( j = 1, 2, JV, ^ = 1, 2, 3), 
the condition of the normalization of the density matrix leads to 



a 



o,o,. ..,o 

1,2 N 



1 



(5) 



Consider an arbitrary spin pair. Without loss of generality we can assume that it includes 
spins 1 and 2. We will investigate the entanglement and the quantum discord for that pair. 
The reduced over spins 3, 4, N density matrix // 1,2 )(t) is 

3 

(6) 



a, 6=o 



where a{]£ 2 (t) = al]& fi 'x'°(t)2 N - 2 . One can verify that Tx{p^} = I. Functions a{^ 2 (t) 
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can be expressed via the correlation functions of the spin system as follows: 

' 4Tr{p(t)(if 1 ® if 2 ® E 3 ® ... g> Ejv)}, 

ei,6 = 1,2,3, 

Tr{p(t)(/f ® if 2 ® £ 3 <8> ... ® Ejv)}, 
6 = 0,6 = l,2,3or ^1 = 1,2,3,6 = 0, 
|Tr{p(t)(/f ® if 2 ® E 3 ® ... ® 



6 = 0,6 = 0. 

One can find from Eq. (jSJ) that a®® = 1/4. The relationships flTj) lead to 



(7) 



^03 ^30 12 21 31 ^13 n 

a 12 = a 12 = a 12 = a 12 = a 12 = a 12 = 0- 



(8) 



Taking into account the permutation symmetry over the numbers of spins 1 an 2, it is 
convenient to introduce the correlation functions 



p = Tr{p(t)(If ®E 2 ®E 3 ®...®E N )} 
= Tr{p(*)(lf ®E x ®Ez® ... <g> E N )}, 
q = TT{p{t){Il g> II g> E 3 ® ••• ® E N )}, 
r = Tr{p(t)(/f ® If ® £ 3 ® ... ® Ejv)}, 
K = Tr{p(t)(Jf ® if (8 S 3 ® ... (8) E N )} 
= Tr {/>(*) (if <g> i 2 2 ® £ 3 ® ... ® E N )}, 
v = Tx{p(t)(If g> I| ® ^3 ® ••• ® £jv)}- 



(9) 

(10) 

(11) 
(12) 

(13) 



Using Eq. (j3J) for the density matrix and executing the unitary transformation U\ = 
exp(i7rJf ) g) _E 2 <g> £7 3 <8> ... <8> £?jv under the sign of the trace (fTBl we obtain that t> = 0. Using 
Eqs. ([7|l- (|13p the reduced density matrix can be written as 



p' 13) (t) 



| — iu f — iu q 



— r\ 



+ iu 



q + r | 



iu 



+ iu q + r 



y g — r 



i 

4 



(14) 



/ 



It is obvious that the density matrix elements p\- (t) = 1,2,3,4) satisfy to the 
relationship pff\t) = p^f^At). Such matrices are called centrosymmetric \x^. One can 
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show that the general expression for the hermitian CS 4x4 matrix is 



/ 



P 



Pi P2 + ip3 PA + iP5 



\ 



P2 ~ip3 \ - Pi Pi Pa ~ ip5 
i 



(15) 



Pa ~ W5 Pi 2~Vi Pi- m 
V Pg Pa + IP5 P2 + IP3 Pi ) 

where pj (i = 1, 7) are real parameters. Because the density matrix must be nonnegative 
defined, all its eigenvalues are greater or equal to zero: 



Ai, 2 = Upe + Pi + 1/2) ± J\{2 Pl +P6-P7- 1/2) 2 + (pa + Pa) 2 + (pa + p 5 ) 2 > 0, 



A 3i4 = |(l/2 -p 6 -p 7 ) ± y^pi -p 6 +p 7 " 1/2) 2 + (p 2 ~ Pa) 2 + (p 3 -p 5 ) 2 > 0. 
These inequalities restrict the range of values for the parameters pi, ...P7. 



(16) 



III. CALCULATION OF THE CORRELATION FUNCTIONS 



To calculate of the correlation functions 



spin 1/2 [20[ and the formula 15] 



-iatll j+ giatrj _ e -iat(2I z -l) j+ 



I3l) we use the commutation relations for 



(17) 



where J + = + As a result, the expressions under the trace sign in Eqs. (J71)-(JT2l) can 
be represented as products of two operators. One of them depends only on spins 1 and 2 
while the other operator depends on spins 3, 4, ... and equals the product of one-spin 
operators. Such approach allows us to rewrite the correlation functions as follows: 

1 



p = - tanh — cos^ 1 (at), 
F 2 2 y 1 



q + r = - tanh 2 — , q — r = — tanh 2 ^ cos^ 2 (2at) 



1 

4 



2 



' - P - .V-2, 



u—- tanh — cos (at) sin (at). 



(18) 
(19) 
(20) 



The correlation functions (fT8|) - (!20|) depend on time periodically and do not decay up to zero 
at t — > oo. This is because the coupling constants for all spin pairs are the same in the 
considered model. 
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IV. ENTANGLEMENT OF SPIN PAIRS IN A NANOPORE 



In our model elements of the CS matrix of Eq. (fT5|) are given by the expressions 
1 p 

Pi = t, P2 = Pi = -, P3 = P5 = -u, Pe = q - r, p 7 = q + r. (21) 

Using Eqs. (1A.9j) one can find that for the system with the density matrix (f!4|) the 
parameters Ai, A2, A3, A4 needed to calculate the concurrence are 



fz fz \ 1 

Ai = — + w, A 2 = — - w, A 3 = - - q - r, A 4 = - - q + r, (22) 

where 

3. 



l + 4g) 2 -16p 2 . - 
, w = vr z + 4u z 



(23) 

Using Eqs. (fl8l) - ( 1201) . ( )22l) . and ( 1231) one can show that Ai is the maximal value in the set 
{Ai, A 2 , A 3 , A 4 }. Then the concurrence (see Eq. (1A.4|) ) is 

C(p (1 ' 2) ) = max jo, 2(Vr 2 + 4m 2 + q) - i | . (24) 

The dependence of the concurrence on the inverse temperatures and dimensionless time 
for different numbers of spins is represented in Fig. [2j It is evident from Fig. |2] that the 
concurrence equals zero at t = 0. Then the concurrence emerges and its evolution has 
periodic character. One can see in Fig. [2] that the concurrence disappears at shorter times 
(inside a period 71/ a) when the number of the spins increases. In particular, one can obtain 
from Eq. (|24|) that the concurrence is absent in the system at iV — > 00. 



V. QUANTUM DISCORD IN SPIN PAIRS IN A NANOPORE 

An investigation of the quantum discord in the system described by the density matrix 
(I14p is a technically complex problem. First we consider the case iV >> 1, at 7^ ttI, 2at 7^ irn 
(I, n are positive integers). Then the correlation functions f 1 1 8 1) - fT20T) are getting simple and 
equal 

p = u = 0, q = r = — tanh 2 — . (25) 

Thus, the density matrix (1141) has the diagonal Bell form for which the discord Q was 
evaluated in [loj]. Using the technique [1Q 1 ], we find that the discord Q for the spin pair in a 
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Figure 2: Concurrence as a function of the dimensionless time r = at and inverse temperatures for 
different numbers of spins N = 6 (a), N = 100 (b). 



nanopore equals 

Q = \ {(1 + 8g) log 2 (l + 8g) + (1 - 8q) log 2 (l - 8q)}-^^ l g 2 (l+4g)-i^ log 2 (l-4g) 

(26) 

At a low temperature T — > (/3 — >■ oo) the asymptotic behavior of the discord is determined 

as 

q^log,--^- (27, 



Notice that for the temperature limit T = 0, (f3 = oo) we have the expression Q = | log 2 | 



0, 3113. Surprisingly, the same analitical result was obtained in [21( for another special mixed 
state. At a high temperature T — > oo (/3 — > 0) the asymptotic behavior of the quantum 
discord is the following 

Q w /3 4 = f-^.V . (28) 

^ 1281n2 1281n2Vfci?Ty v ; 

The dependence of the quantum discord on the temperature is represented in Fig. [3j One 
can see from Fig. |3] that quantum correlations are relatively large only at milli-kelvin 
temperatures for the Larmour frequency Uq = 2tt ■ 500 ■ 10 6 s _1 . 
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Figure 3: The dependences of the quantum and geometric discord for a spin pair in a nanopore on 
the temperature at iV >> 1, at ^ ttI, 2at ^ ttti (l,n are positive integers). The Larmour frequency 
u = 2vr • 500 • 10 6 a -1 . 

The considered model allows us to investigate some peculiarities of the time evolution of 
the discord. At the time moments U = (1 + 2Z)7r/(2a), Z = 0,1,2... the correlation functions 
are 

P = u = 0, q = I tanh 2 t [1 + (_l)^] , r = I tanh 2 | [l - (-l)^ 2 ] . (29) 

Essentially that one of the correlation functions, q or r, equals zero depending on the parity 
of the number of spin-carrying molecules in the nanopore. Here the structure of the density 
matrix coincides with so-called X-matrix for which the method of the calculation of the 
discord was worked out in Q. Tho method Q allows us to eouclude that the quautum 
discord equals zero at the time moments ti = (1 + 21)tt/ (2a) (I = 0, 1, 2, ...). Such is the case 
because one of the correlation functions equals zero in these time moments. Thus we have 
"flickering" quantum discord which equals zero periodically. 



VI. GEOMETRIC DISCORD FOR THE 4 x 4 CS DENSITY MATRIX 



The geometric discord 



22j simplifies calculations of the quantum discord. In particular, 



it allows us to find the discord for the CS matrix (Jl5j) . Using the Pauli matrices a) 
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21 j (j = 1, 2; i = x, z) we can rewrite the density matrix (|T5|) in the Bloch form 



i + £ T ij( r i <g> ^' + ELi ^ ® 1 + Eti y»i ® ^ 



= I [1 + (A Pl - l)a z ®a z + 2(p 6 + p 7 )o x ® a x + 2(p 7 - p 6 )aV g> <r»- 

-4p 3 a z ®a y - Ap 5 a y ® a z + Ap A a x ® 1 + Ap 2 l ® a x ) , 
where T is the correlation matrix, equal to 

f 2(p 6 +p r ) ^ 

T= 2(p 7 -p 6 ) -4p 5 , 

v -Ap 3 4pi - ly 



and x, 1/ are the Bloch vectors 



5 T = (4p 4 , o, o) , y T = (4p2, o, o) . 



The geometric discord Q 9 for the two-qubit system is given by [22] 

= g (ll ^ II ll-^ll — krnax^j i 

where the spheric (Hilbert-Schmidt) norm \\A\\ is 



Tr (AA+) 



and k max is the maximal eigenvalue of the matrix 

— >— >^ rp 

K=xx +TT . 

The geometric discord can be expressed through ki,k 2 , k 3 as 

1 

Qg = 2 + ^2 + ^3 — ^mi) • 

Simple calculations yield the eigenvalues ki,k 2 , k 3 of the matrix K 

h = 16pl + A{p 6 +p 7 ) 2 , 
k 2 , 3 = 8pl + 2(p 7 - p 6 ) 2 + 8p| + |(4 Pl - 1) 2 ± 
± {[8p 2 + 2 (p 7 - p 6 ) 2 - 8p 2 - i(4 Pl - l) 2 ] 2 + 16[p 5 (4 Pl - 1) + 2p 3 (p 7 - Pe)] 2 }^ 
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In case of a nanopore we can obtain using (I2TI) the expressions 

h = 4(jo 2 + 4g 2 ), 



, (38) 

k 2 , 3 = 8 [2u 2 + r 2 ± ry/r 2 + 4u 2 \ . 

As a result, the geometric discord has the form 

Q, = ltanh 4 ^. (39) 

In the limit T — y 0, the geometric discord Q g = 1/8 = 0.125. At a high temperature 
T — > oo, (f3 — > 0) the asymptotic behavior of the geometric discord is the following 



1 _ 4 1 / fojQ 



i 



Q q « — /3 4 = — ( f^- ] . (40) 
^ 9 128 128 \k B T 1 



VII. CONCLUSION 



We investigated the entanglement and quantum discord in a nanopore filled with a gas 
of spin-carrying molecules (atoms). Even when the entanglement is absent the quantum 
discord can reach large values in the considered system. It reveals significant quantum 
correlations in the system. They exhibit "flickering" character and equal zero periodically 
in the process of the time evolution of the system. 

The density matrix in the considered model belongs to the so-called CS matrices. We 
evaluated the analytic expression for the concurrence of a general two-qubit CS density 
matrix. We studied the geometric discord for systems with CS density matrices and obtained 
an analytic expression for the geometric discord in the considered model. 

The model under question confirms once more that entanglement describes only a part 



of quantum correlations while the discord is a measure of total quantum correlations [l]. 

The authors thank A.I. Zenchuk for useful discussions. The work is supported by the 
Program of the Presidium of RAS No. 8 "Development of methods of obtaining chemical 
compounds and creation of new materials". 



Appendix. Entanglement of the CS 4x4 density matrix 



To calculate the entanglement for the general two-qubit density matrix, one should find 
the "spin-flip transformed" density matrix which is 18 j 



p=K®a 2 >*K®af), 



(Al) 
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where p* is the complex conjugate matrix to p in the standard basis {| 00) , |01) , |10) , |11)}, 
and a 3 y (j — 1, 2) is the Pauli matrix. Then it is necessary to calculate the concurrence C(p) 
which determines unambiguously the entanglement E(p) 18 1 

l + y/l-C*(p) 



E{p) = H 

where H(x) is the Shannon function 23j 



H(x) = —x log 2 x — (1 — x) log 2 (l — x). 
The concurrence is given as [3] 

C{p) = max {0, 2A - Ai - A 2 - A 3 - A 4 } 



(A.2) 



(A.3) 



(A.4) 



A = max {Ai, A 2 , A 3 , A 4 } , 

where Ai, A 2 , A 3 , A 4 are the square roots of the eigenvalues of the matrix product pp. 

One can show that the matrix p is a CS matrix if the density matrix p is a CS one. It 
means that the matrix pp is also a CS one. Under the orthogonal transformation 

/l 1 \ 



V2 



S 1 



(A.5) 



Oil 
1-10 
\1 -1/ 

any 4x4 CS matrix takes the block- diagonal form consisting of two 2x2 subblocks. As a 
result, the square roots of the eigenvalues of the matrix product pp can be written as follows 

Ai, 2 = \ { (pi + p 6 ) 2 - 2(p 2 + p 4 ) 2 + 2(p 3 + p 5 ) 2 + (§ - pi + vif 



(A.6) 



±yJ[(2 Pl + p 6 - \ - p 7 ) 2 + 4(p 3 + p 5 ) 2 ] [(| + p 6 + p 7 ) 2 - 4(p 2 + p 4 ) 2 ] 
A 3 ,4 = 75 { (Pi - Pe) 2 - 2(p 2 - p 4 ) 2 + 2(p 3 - p 5 ) 2 + (| - pi - p 7 ) 2 

±^/[(2 Pl - p 6 - \ + p 7 ) 2 + 4(p 3 - p 5 ) 2 ] [(i - p 6 - p 7 ) 2 - 4(p 2 - p 4 ) 2 ] } 
Using identities 

2 [(Pi + Pe) 2 + (| - Pi + Pi) 2 } = (2pi + p 6 - | - p 7 ) 2 + (| + Pe + p 7 )' 
2 [(Pi - Pe) 2 + (| - Pi ~ Pi) 2 } = (2pi - Pe - | + p 7 ) 2 + (| - Pe - Pr) 5 



(A.7) 



(A.8) 
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we can rewrite finally Ai, A2, A3, A4 as 



Al = \ { yJ(2 Pl +P6-|- P7 ) 2 + 4( P3 + P 5 ) 2 + ^(i+P6+P7) 2 -4(p 2 +p 4 ) 2 } , 



A 2 — I 



(2pi +p 6 -|- p 7 ) 2 + 4(p 3 + P5 y -Ja+p 6 + P7 y - 4{ P2 + P4 y 



(A.9) 



A 3 



\ { yJ(2 Pl -p6-i+ p 7 ) 2 + 4(p 3 - Ps) 2 + - Pe - Pr) 2 - 4(p 2 - p 4 ) 2 } 



A 4 = i ^(2^ - p 6 - i + P7 y + 4(p 3 - P5) 2 - - Pe - Pr) 2 - 4(p 2 - p 4 ) 2 • 
These relationships together with Eqs. (IA.2j) - flA.4j) give an analytical formula only through 
the square radicals for calculating the entanglement of formation for the arbitrary two-qubit 
density matrix (I15p. 
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